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. A random-field Ising model that is capable of exhibiting a rich variety of multi- 

I , critical phenomena, as well as a smearing of such behavior, is investigated. The 

^ ' model consists of an infinite-range-interaction Ising ferromagnet in the presence of a 

o ■ 

, triple-Gaussian random magnetic field, which is defined as a superposition of three 
Gaussian distributions with the same width cr, centered at = and H = HJq, 

> 

10 ' with probabilities p and (1 — p)/2, respectively. Such a distribution is very general 

00 ■ 

• and recovers as limiting cases, the trimodal, bimodal, and Gaussian probability dis- 

CN 

. tributions. In particular, the special case of the random-field Ising model in the pres- 

00 ! ence of a trimodal probability distribution (limit a ^ 0) is able to present a rather 

o 

nontrivial multicritical behavior. It is argued that the triple-Gaussian probability 



^ ■ distribution is appropriate for a physical description of some diluted antiferromag- 

H : 

nets in the presence of a uniform external field, for which the corresponding physical 
realization consists of an Ising ferromagnet under random fields whose distribution 
appears to be well-represented in terms of a superposition of two parts, namely, a 
trimodal and a continuous contribution. The model is investigated by means of the 
replica method, and phase diagrams are obtained within the replica-symmetric so- 
lution, which is known to be stable for the present system. A rich variety of phase 
diagrams is presented, with a single or two distinct ferromagnetic phases, continuous 
and first-order transition lines, tricritical, fourth-order, critical end points, and many 
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other interesting multicritical phenomena. Additionahy, the present model carries 
the possibihty of destructing such multicritical phenomena due to an increase in 
the randomness, i.e., increasing a, which represents a very common feature in real 
systems. 
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I. INTRODUCTION 



The Random Field Ising Model (RFIM) became nowadays one of the most studied 
problems in the area of disordered magnetic systems hj, ^1. From the theoretical point 

n n 

of view [3|, its simple definition ^4!], together with the richness of physical properties that 
emerge from its study, represent two main motivations for the investigation of this model. 
On the other hand, a considerable experimental interest js] arised after the identification 
of the_^RFIM with diluted antiferromagnets in the presence of a uniform magnetic field 
6 

and Fe,Mgi_,Cl2 |5|,l9|. 

In what concerns the equilibrium phase diagrams of the RFIM, the effects of differ- 
ent probability-distribution functions (PDFs) for the random fields have attracted the 



,171,1811; since then, two of the most investigated systems are the compounds FexZni_xF2 
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authors, see, e.g., Refs. [10 
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26| . among others. At the mean-field level, the Gaussian PDF yields a 



continuous ferromagnetic-paramagnetic boundary [lOj, whereas discrete PDFs may lead 
to elaborate phase diagrams, characterized by a finite-temperature tricritical point fol- 

Q, Q , or even 
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lowed by a first-order phase transition at low temperatures 

fourth-order and critical end points [l^, isl. For the short-range- interaction RFIM, the 
existence of a low-temperature first-order phase transition remains a very controversial 
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From the experimental point of view, one is certainly concerned with physical realiza- 
tions of the RFIM, for which the most commonly known are diluted antiferromagnets in 



the presence of a uniform magnetic field 



5|]. Hence, in the identification of RFIMs with 



such real systems , certain classes of distributions for the random fields, in the 

corresponding RFIM, are more appropriate for a description of diluted antiferromagnets 
in the presence of a uniform magnetic field. In the later systems one has local variations 
in the sum of exchange couplings that connect a given site to other sites, leading to lo- 
cal variations of the two-sublattice site magnetizations, and as a consequence, one may 
have local magnetizations that vary in both sign and magnitude. In the identifications 
of the RFIM with diluted antiferromagnets 0, Q, , the effective random field at a given 
site is expressed always in terms of quantities that vary in both sign and magnitude: 
(i) the local magnetization jol; (ii) two contributions, namely, a first one that assumes 
only three discrete values, related to the dilution of the system and the uniform external 
field, and a second one that is proportional to the local magnetization [3]; (iii) the sum 
of the exchange couplings associated to this site [8]. Therefore, for a proper description 
of diluted antiferromagnets in the presence of a uniform field, the corresponding RFIMs 
should always be considered in terms of continuous PDFs for the fields. 

The compound FexMgi_xCl2 presents an Ising spin-glass behavior for x < 0.55, and is 
considered as a typical RFIM for higher magnetic concentrations. In the RFIM regime, it 
shows some curious behavior and is considered as a candidate for exhibiting multicritical 
phenomena . As an example of this type of effect, one finds a first-order transition 



turning into a continuous one due to a change in the random fields jsl, 128 1; the con- 
centration at which the first-order transition disappears is estimated to be a; = 0.6. The 
crossover from first-order to continuous phase transitions has been investigated through 



different theoretical approaches 



13, 



28|. One possible mechanism used to find such a 
crossover, or even to suppress the first-order transition completely, consists in introduc- 
ing an additional kind of randomness in the system, e.g., bond randomness 16|, ll7|. By 
considering randomness in the field only, this crossover has been also analyzed through 
zero-temperature studies, either within mean- field theory 29|], or numerical simulations 



on a three-dimensional lattice 



28 



2^ ■ Recently a RFIM has been proposed 26j for which 
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the finite-temperature tricritical point, together with the first-order line, may disappear 
due to an increase in the field randomness, similarly to what happens with the first-order 
phase transition in the compound FexMgi_xCl2. 

However, it is possible that the diluted antiferromagnet FexMgi_xCl2, or some other 
similar compound, may present an even more complicated critical behavior, not yet ver- 
ified experimentally, to our knowledge. According to the analysis of Ref. [7|, one of its 
contributions for the random fields assumes only the values 0, ±\^H (H represents the 

nn 

external uniform magnetic field); motivated by this result, a RFIM was proposed [141. Il5l| 
with the random fields described in terms of a trimodal distribution. Such a RFIM, stud- 
ied within a mean-field approach through a model defined in the limit of infinite-range 
interactions, yielded a rich critical behavior with the occurence of first-order phase tran- 
sitions, tricritical and high-order critical points 30l at finite temperatures, and even the 

nn 

possibility of two distinct ferromagnetic phases at low temperatures [U, 115||- These inves- 
tigations suggest that FexMgi_xCl2 may exhibit a rich critical behavior that goes beyond 
the disappearance of a first-order phase transition observed around x = 0.6. However, 
taking into account the above criteria for the identification of the RFIM with diluted 
antiferromagnets, one notices that the trimodal distribution does not represent an ap- 
propriate choice from the physical point of view. In fact, in the analysis of Ref. [t], a 
second, continuous contribution for the random fields should be taken into account as 
well. Therefore, for an adequate theoretical description of FexMgi_xCl2, one should con- 
sider a RFIM defined in terms of a finite-width three-peaked distribution, which could 
exhibit multicritical behavior, and in addition to that, a crossover from such behavior to 
continuous phase transitions due to an increase in the field randomness. 

For that purpose, herein we introduce a RFIM that consists of an infinite-range- 
interaction Ising ferromagnet in the presence of a triple- Gaussian random magnetic field. 
This distribution is defined as a superposition of three Gaussian distributions with the 
same width a, centered at if = and H = ±Hq, with probabilities p and (1 — p)/2, 
respectively. T he p articular case of the RFIM in the presence of a trimodal PDF, studied 
previously [l4, isf, is recovered in the limit a 0. We show that the present model is 
capable of exhibiting a rich multicritical behavior, with phase diagrams displaying one, or 
even two distinct ferromagnetic phases, continuous and first-order critical frontiers, tri- 
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critical points and high-order critical points, at both finite and zero temperatures. Such 
a variety of phase diagrams may be obtained from this model by varying the parameters 
of the corresponding PDF, e.g., p and a: (i) Variations in p {a fixed) lead to qualitatively 
distinct multicritical phenomena; (ii) Increasing a {p fixed) yields a smearing of an specific 
type of critical behavior, destroying a possible multicritical phenomena due to an increase 
in the randomness, representing a very common feature in real systems. In the next sec- 
tion we define the model, find its free energy and equation for the magnetization, by using 
the replica approach. In section III we present several phase diagrams, characterized by 
the critical behavior mentioned above. Finally, in section IV we present our conclusions. 



II. THE MODEL 



The infinite-range-interaction Ising model in the presence of an external random mag- 
netic field is defined in terms of the Hamiltonian 

where the sum j) ^^ns over all distinct pairs of spins Si — ±1 (i = 1,2, ...,A^). The 

random fields {Hi\ are quenched variables and obey the PDF, 



{1-p) 



1 



1/2 



+ P 



27r(72 



27r(72 

1/2 

exp 



exp 



[H^ - Ho) 

2(72 



+ exp 



[H, + Ho) 

2(72 



H? 



2(72 



(2.2) 



which consists of a superposition of three independent Gaussian distributions with the 
same width a, centered at Hi = and Hi = ±Ho, with probabilities p and (1 — p)/2, 
respectively (to be called hereafter, a triple Gaussian PDF). This distribution is very 
general, depending on three parameters, namely, p, a, and Hq, and contains as particular 
cases several well-known distributions of the literature, namely the trimodal and bimodal 
distributions, the double Gaussian, as well as simple Gaussian distributions. The model 
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defined above is expected to be appropriate for a pliysical description of diluted antifer- 
romagnets tliat in the presence of an external magnetic field may exhibit multicritical 
behavior; one candidate for this purpose is FexMgi_xCl2 0]- 

From the free energy F{{Hi}), associated with a given reahzation of site fields {Hi}, 
one may calculate the quenched average, [F{{Hi})]H, 

mmu = I x{mp{Hi)]Fm}) . (2.3) 

The standard procedure for carrying the average above is by making use of the replica 
method {2, 31], leading to the free energy per spin, 



f3f = hm h\nZi{H,})]H = lim lim -^([Z^^ - 1) , (2.4) 



where is the partition function of n copies of the original system defined in Eq. (12.11) 
and P = l/{kT). One gets that, 

/?/ = lim — min g{m") , (2-5) 

n^O n 

with 

ai^"") = ^E(^")'-^^l^Tr.^xp(7^+)-ii^lnTr„exp(7^;J 

a 

- plnTr„exp(HS)) , (2.6) 
K'iF = /^-^E^^^^ + Z^^lE^J ±PHoJ2s\ (2.7) 

a \ a / a 

(2.8) 



a \ a / 



In the equations above, a represents a replica label {a = 1,2, ...,n) and Tr^ stands for a 
trace over the spin variables of each replica. The extrema of the functional g{m°') leads 
to the equation for the magnetization of replica a, 
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= < 5" >+ + ^^r-^ < 5" >_ + p < 5° >o , (2.9) 



where < >± and < >o denote thermal averages with respect to the "effective Hamiltoni- 
ans" H^ff and Tig^, in Eqs. (12.71) and (12. Sp . respectively. 



If one assumes the replica-symmetry ansatz [2, l3l|, i-e., m'^ = m (V a), the free energy 



per spin of Eqs. (I2.5p - (I2.8I) and the equilibrium condition, Eq. (12. 9p . become 

/ = ^m^- j D2ln(2cosh$+) - ^^^^^ j ln(2 cosh $~) 

- ly D2ln(2cosh$(°)) , (2.10) 

m = ilzi^ / Dz tanh<i>+ + 2 / ^^^a^h'^^+^P y ^^tanh$(°) , (2.11) 



where 



/I r+oo 
Dz (..) = ^= / c/ze-''/2 (..) . 

(2.12) 

i should be mentioned that the instability associated with the replica-symmetric solution 
32] at low temperatures is usually related to parameters characterized by two replica 



indices, like in the spin-glass problem [2, l3ll- In the present system the order parameter 



m" depends on a single replica index, and so such an instability does not occur. 

In the next section we will make use of Eqs. (I2.10p and (12.111) in order to get several 
phase diagrams for the present model. Although most critical frontiers will be achieved 
numerically, some analytical results may be obtained at zero temperature. At T = 0, the 
free energy and magnetization become, respectively. 
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(2.13) 



III. PHASE DIAGRAMS 

Usually, in the RFIM one has two phases, namely, the Ferromagnetic (F) and Param- 
agnetic (P) ones. However, for the PDF of Eq. (12. 2p . there is a possibility of three distinct 
phases: apart from the Paramagnetic (m = 0), one has two Ferromagnetic phases, Fi and 
F2 , characterized by different magnetizations (herein we will consider always Fi as the 
ordered phase with higher magnetization, i.e., mi > m2 > 0). Close to a continuous 
transition between an ordered and the disordered phases, m is small, so that one can 
expand Eq. (12.111) in powers of m, 

m = Aim + A^m^ + A^m^ + A-juiJ + O(m^) , (3.1) 
where the coefficients are given by 



Ai = (3J{l-(l-p)X,-p\f^} , (3.2) 

As = -^{l-4[(l-p)Ai+rf)]+3[(l-p)A2+rf)]}, (3.3) 

A, = i^{2-17[(l-p)Ai+rf)] + 30[(l-p)A2+pAf]-15[(l-j9)A3 
15 

+ rf)]} (3.4) 



A7 = ^{17- 248[(l-j9)Ai+pAH + 756[(l-p)A2 + pAf]- 840[(l-p)A3 
+ pAf]+315[(l-p)A4+pAf]}, (3.5) 
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with 



'2tx J 



/ + 00 
dz e'^'/2 ^^^^2k p^jj^ _^ ^^-^ 
-oo 



oo 
+00 



1 p+co 

\f = -j= / dze-'"/Hanh''\(3az) . (3.6) 

In order to find the continuous critical frontier one sets Ai = 1, provided that A3 < 0. If 
a first-order critical frontier also occurs, the continuous line ends when ^3 = 0; in such 
cases, the continuous and first-order critical frontiers meet at a tricritical point, whose 
coordinates may be obtained by solving the equations Ai = 1 and ^3 = numerically, 
provided that < 0. In the present problem there is also a possibility of a fourth-order 
critical point, which is obtained from the conditions Ai = 1, A3 = A^ = and Aj < 0. 
In addition, when the two distinct ferromagnetic phases are present, two other critical 
points may also appear (herein we follow the classification due to Griffithis jsol): (i) the 
ordered critical point, which corresponds to an isolated critical point inside the ordered 
region, terminating a first-order line that separates phases Fi and F2 ; (ii) the critical end 
point, where all three phases coexist, corresponding to the intersection of a continuous 
line that separates the paramagnetic from one of the ferromagnetic phases with a first- 
order line separating the paramagnetic and the other ferromagnetic phase. The location 
of the critical points defined in (i) and (ii), as well as the first-order critical frontiers, were 
determined by a numerical analysis of the free-energy minima, e.g., two equal minima for 
the free energy, characterized by two different values of magnetization, mi > 1712 > 0, 
yields a point of the first-order critical frontier separating phases Fi and F2 . 

Next, we show several phase diagrams of this model for both finite and zero tempera- 
tures. Besides the notation defined above for labeling the phases, in these phase diagrams 
we shall use distinct symbols and representations for the critical points and frontiers, as 
described below. 

• Continuous critical frontier: continuous line; 

• First-order critical frontier: dotted line; 

• Tricritical point: located by a black circle; 
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• Fourth-order point: located by an empty square; 

• Ordered critical point: located by a black asterisk; 

• Critical end point: located by a black triangle. 

A. Finite- Temperature Critical Frontiers 

We now present finite-temperature phase diagrams; these phase diagrams will be ex- 
hibited in the plane of dimensionless variables, {kT)/J versus Hq/J, for typical values of 
p and a/ J. For completeness, in each figure we also show the forms of the corresponding 
PDF of Eq. (12.21) . with the parameters used in the phase diagrams. 

In Fig. [1] we present some qualitatively distintic phase diagrams of the model, for the 
fixed value p = 0.2 and typical values of a/ J. In Fig. [^a) we represent the PDF of 
Eq. (12. 2p for p = 0.2, (Hq/J) = 0.56, and the widths a/ J used to obtain the phase 
diagrams of Figs, [m^c)-(f). The particular case (cr/J) = 0.1, used in Fig. [^b), yields a 
three-narrow-peaked distribution and is not represented in Fig. [T](a) for a better visualiza- 
tion of the remaining cases; in addition to that, the choice [Hq/ J) = 0.56 corresponds to a 
region in the phase diagrams where essential changes occur in the criticality of the system 
at low temperatures. One notices that the critical frontier separating the paramagnetic 
(P) and ferromagnetic (Fi and F2) phases exhibits significant changes with increasing 
values of a/ J. For (cr/J) = 0.1 [cf. Fig. [11(b)], one has a phase diagram that is quali- 



tatively similar to the one obtained in the case of the trimodal distribution [15|, with a 
tricritical point (black circle) and a critical end point (black triangle) at a lower, but finite 
temperature, where the critical frontier separating phases P and F2 meets the first-order 
line. However, for increasing values of a/ J these critical points move towards low tem- 
peratures, in such a way that for (cr/J) = 0.16 one observes the collapse of the critical 
end point with the zero-temperature axis, i.e., the ferromagnetic phase F2 occurs only for 
T = 0, as shown in Fig. [11(c). Therefore, for p = 0.2 the value (cr/J) = 0.16 represents a 
threshold for the existence of phase F2. One should notice from Fig.[Tl(a), that for p = 0.2, 
the existence of the phase F2 is associated with a PDF for the fields characterized by a 
three-peaked shape. For (a/ J) = 0.3 [cf. FigJHd)], the phase diagram is qualitatively 



similar to the one obtained for the bimodal 111 



or continuous two-peaked distributions 




FIG. 1: (a) The probability distribution of Eq. (|2.2p . for p = 0.2 and {Hq/ J) = 0.56 is represented 
for several values of its width, (b)-(f) Phase diagrams showing the critical frontiers separating 
the paramagnetic (P) and ferromagnetic (Fi and F2) phases of the infinite-range-interaction 
ferromagnet in the presence of a triple Gaussian random field, for typical values of cr/J and 
p = 0.2. Critical frontiers and critical points are as described in the text. All quantities are 
scaled in units of J. 



26l ] , where one finds a tricritical point at finite temperatures and a first-order transition at 
lower temperatures. We have found analitically through a zero-temperature analysis that 
will be discussed later, the value (cr/J) = 0.4551 for which the tricritical point reaches the 
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FIG. 2: (a) The probability distribution of Eq. ([22D, for p = 0.308561 and {Hq/J) = 0.65, is 
represented for several values of its width, (b)-(d) Phase diagrams showing the critical frontiers 
separating the paramagnetic (P) and ferromagnetic (Fi and F2) phases of the infinite-range- 
interaction ferromagnet in the presence of a triple Gaussian random field, for typical values of 
a/ J and p = 0.308561. Critical frontiers and critical points are as described in the text. All 
quantities are scaled in units of J. 

zero-temperature axis, signaling a complete destruction of the first-order phase transition 
in the case p = 0.2, as exhibited in Fig. [11(e). In this case, the associated PDF presents a 
single flat maximum in agreement with Refs. ll|, ll2|, ll3| . For even higher values of cr/J, 
as shown in Fig. [T](f) for the case (cr/J) = 0.6, the PDF presents a single peak, and the 
frontier ferromagnetic-paramagnetic is completely continuous, qualitatively similar to the 



case of a single Gaussian PDF 10 1 



There is an special value of p, to be denoted herein p*, which represents an upper limit 
for the existence of tricritical points along the paramagnetic border, characterized by a 
fourth-order point at zero temperature. The value p* = 0.308561, as well as its associated 
width, {a* /J) = 0.369294, will be determined later, within a zero-temperature analysis. 
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In Fig. [2] we show phase diagrams for p = 0.308561 and three different values of a/ J. In 
Fig. [2](a) we represent the PDF of Eq. (Q for p = 0.308561, {Hq/J) = 0.65, and the 
widths a/ J used to obtain the phase diagrams of Figs. [2]^b)-(d); the choice {Hq/ J) = 0.65 
has to do with a region in the phase diagrams where interesting critical phenomena occur 
at low temperatures. In Fig. EJ^b) we present the phase diagram for (cr/J) = 0.1, where 
one sees that the border of the paramagnetic phase is completely continuous; for lower 
(higher) fields this critical frontier separates phases P and Fi (F2). However, there is a 
curious first-order line separating phases Fi and F2 that terminates in an ordered critical 
point (black asterisk), above which one can pass smoothly from one of these phases to the 
other. By slightly increasing the values of a/ J one observes that the ordered critical point 
disappears, leading to the emergence of both tricritical and critical end points at finite 
temperatures, yielding a phase diagram that is qualitatively similar to the one shown in 
Fig. [U^b); this type of phase diagram occurs typically in the range 0.1 < (cr/J) < 0.3. 
However, as shown in Fig. [2]^c) for (cr/J) = 0.305, one has the collapse of the critical end 
point with the zero-temperature axis, with the phase F2 occurring at zero temperature 
only, analogous to the one that appears in Fig. [^c), but now for a higher value of Hq/J. 
Comparing Figs. [H^c) and M^c), one sees that in some cases, essentially similar phase 
diagrams may be obtained by increasing both p and cr/J; however, in the case of higher 
values for these parameters, the extension of the first-order transition line gets reduced. It 
should be mentioned that, in both figures, the existence of the phase F2 is associated with 
a PDF for the fields characterized by a three-peaked shape. In Fig. [2](d) we display the 
phase diagram for (cr/J) = 0.369294, for which the first-order line is totally destroyed. It 
should be emphasized that this effect appears also for other values of p through a tricritical 
point at T = [cf. Fig. W^e)]; however, the phase diagram of Fig. [2](d) is very special, 
in the sense that the collapse of the first-order frontier with the zero-temperature axis 
occurs by means of a fourth-order point. For higher values of cr/J, the critical frontier 
separating phases P and Fi is completely continuous, and one has essentially the same 
phase diagram shown in Fig. [l|f). 

Additional phase diagrams are shown in Fig. [3l for the case p = 0.5. In Fig. [3t^a) 
we represent the PDF of Eq. (12. 2p for such a value of p, {Hq/J) = 0.74, and the widths 
cr/J used to obtain the phase diagrams of Figs. [3](b)-(d); the choice {Hq/J) = 0.74 
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FIG. 3: (a) The probability distribution of Eq. ([22]), for p = 0.5 and {Hq/J) = 0.74, is rep- 
resented for several values of its width, (b)-(d) Phase diagrams showing the critical frontiers 
separating the paramagnetic (P) and ferromagnetic (Fi and F2) phases of the infinite-range- 
interaction ferromagnet in the presence of a triple Gaussian random field, for typical values of 
a/ J and p = 0.5. Critical frontiers and critical points are as described in the text. All quantities 
are scaled in units of J . 

corresponds to a region of the phase diagram where an ordered critical point appears at 
low temperatures. In Fig. [3t^b) we present the phase diagram for the same a / J value of 
Figs. [D^b) and[2](b), i.e., (cr/J) = 0.1; these figures suggest that the existence of phase 
F2 is associated with three-peaked distributions. However, for sufficiently small p, this 
phase appears together with a critical end point, as shown in Fig. [T](b), whereas for larger 
p, the border of the paramagnetic phase is completely continuous, and phases Fi and 
F2 are separated by a first-order critical frontier that terminates in an ordered critical 
point. In Fig. [3](b) one can go smoothly from one of these two ferromagnetic phases 
to the other, through a thermodynamic path connecting these phases above the ordered 
critical point. In Fig. [3](c) one has the collapse of the ordered critical point with the 
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FIG. 4: Projections of the fourth-order-point hne on two different planes: (a) The plane {kT)/J 
versus p, showing the critical temperatures associated with fourth-order points that exist in 
the interval 8/33 < p < p*; (b) The plane a/ J versus p. The empty squares are computed 
fourth-order points, whereas the lines (guides to the eye) represent lines of fourth-order points. 

zero-temperature axis, which was estimated numerically to occur for (cr/J) = 0.20, with 
the phase F2 appearing at zero temperature. It is important to mention that the phase 
diagram exhibited in Fig. [3](c) is qualitatively distinct from those of Figs. [U^c) and[2]^c), 
in the sense that the first one is characterized by a border of the paramagnetic phase 
that is completely continuous and at zero temperature, phases Fi and F2 are separated 
by an ordered critical point. For the value (cr/J) = 0.4 shown in Fig. [21(d), one gets a 
continuous critical frontier separating phases P and Fi. However, there is a basic difference 
between the phase diagrams displayed in Figs. [T|^f) and[3]^d): as it will be shown in the 
following zero-temperature analysis, in the later case there is no zero-temperature point 
separating these two phases, i.e., the phase Fi exists for all values of Hq/J. In fact, for any 
(cr/J) < 0.4 the border of the paramagnetic phase never touches the zero-temperature 
axis, as suggested in Figs. [31^b)-(d); for higher values of cr/J, this critical line meets the 
Ho/ J axis, and one has a phase diagram qualitatively similar to the one of Fig. [11(f). 

There are some important threshold values associated with the existence of the above- 
mentioned critical points, as described next. 

(i) We found numerically that for p ^ p*, critical end points do not occur, either in 
finite or zero temperatures. 

(ii) We verified, also numerically, that for p > 0.93, ordered critical points cease to 
exist; as a consequence, this represents a threshold for the existence of phase F2, which 
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does not appear above this value (even at zero temperature). 

(iii) In the zero-temperature analysis (to be carried below), we find analytically that 
for p* = 0.308561 one gets a fourth-order point at zero temperature; this value depicts an 
upper bound for the existence of tricritical points. 

(iv) Fourth-order points usually delimitate the existence of tricritical points and are 
sometimes considered in the literature as "vestigial" tricritical points [l5l]. In the case 
of the power expansion in Eq. fl3.ip . they exist for finite temperatures as well, and are 
determined by the conditions Ai = 1, = = and A^ < 0, which define a line 
in the four-dimensional space [{kT)/ J,p, Hq/ J,a/ J]. In Fig. H] we present projections of 
the fourth-order-point line with the planes {kT)/J versus p [Fig. 111(a)] and a/ J versus p 
[Fig. 111(b)]. These projections interpolate between the fourth-order points occurring for 
p = 8/33 (cr = 0, i.e., trimodal distribution |l5^]) and the zero-temperature threshold value 
for the triple Gaussian PDF, p* = 0.308561, to be determined below. 

From the analysis above one concludes that it is possible to obtain qualitatively similar 
phase diagrams for different pairs of parameters {p,cr/J). Essentially, these phase dia- 
grams are defined by the presence of the different types of critical points that may appear 
in this model. The existence of tricritical points was already discussed in items (iii) and 
(iv) above, whereas the regions in the plane a/ J versus p associated with ordered and 
critical end points are exhibited in Fig. (Sj^a). Along the axis (a/J) = our limits for the 
existence of critical end points (0 < p < 0.24) and ordered critical points (0.27 < p < 0.93) 
are in agreement with those estimated in Ref. [isl; in between these two regimes, a small 
region occurs, which is very subtle, from the numerical point of view, where two critical 
end points and one ordered critical point show up. This intermediate region gets reduced 
for increasing values of a/J and is delimited by the two dashed lines shown in Fig. [51(a). 
A typical phase diagram inside this region is shown in Fig. [51(c), which is characterized by 
a critical frontier separating phases P and Fi that presents a continuous piece, followed 
by a first-order line at lower temperatures, with no tricritical point. It is important to 
stress that the type of phase diagram of Fig. [51(c) occurs in the particular case (a/J) = 
for a range of values of p right above p = 8/33 ~ 0.24, which represents the probability at 
which a fourth-order point occurs at finite temperatures |l5|. We have verified that a sim- 
ilar effect occurs for finite values of a/ J, in the sense that the narrow intermediate region 
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FIG. 5: (a) Regions, in the plane a /J versus associated with quahtatively distinct phase 
diagrams for the present model; the empty squares represent computed points, whereas the lines 
are just guides to the eye; the arrow indicates the threshold value p = p*. In regions I and 
II one has phase diagrams like those shown in (b) and (d), respectively. The two dashed lines 
that appear in (a) define a narrow intermediate region, exhibiting a phase diagram characterized 
by two critical end points and one ordered critical point along the paramagnetic border, like 
shown in (c). The type of phase diagram in the intermediate region appears when one varies 
the parameters of the distribution of Eq. ()2.2[) in such a way to go from region I to region II [as 
shown typically in the sequence (b), (c) and (d), for {a/ J) = 0.1], and vice versa. 

of Fig. [5](a) corresponds to a region of points to the right of the projection of the fourth- 
order line in the plane a/ J versus p of Fig.|4](b). This intermediate region occurs whenever 
the fourth order point appears for finite temperatures and it becomes essentially undis- 
cernible, from the computational point of view, as one approaches the threshold p = p*, 
corresponding to the fourth-order point at zero temmperature; we verified numerically 
that this region disappears completely at the threshold p = p* [this aspect is reinforced 
in Fig. [5](d) where we present a phase diagram for p slightly larger than p = p*]. The fact 



18 



that the intermediate region gets reduced for increasing values of a is expected, since one 
of the effects introduced by the parameter a is to destroy several types of multicritical 
behavior. In Fig. [5](a), region I is associated with phase diagrams that present a single 
critical end point, like the one of Fig. [5t^b), whereas region II is associated with phase 
diagrams presenting a single ordered critical point, as shown in Fig.[5](d); one may go from 
region I to region II through the sequence of phase diagrams presented in Figs. E](b)-(d). 
The ranges of p values associated with these two regions diminish for increasing values 
of 0"/J, as shown in Fig. [5](a); the computed points along the full lines (empty squares) 
represent upper limits (in a/ J) for regions I and II, in which cases the corresponding 
critical points appear at zero temperature. As typical examples, one has the ranges for 
the appearence of these points, a/ J < 0.08 {p = 0.1) and a / J < 0.16 {p = 0.2), in region 
I, whereas cr/J < 0.2 {p = 0.5) and cr/J < 0.127 {p ~ 0.6667), in region II. Therefore, 
for a given pair of parameters {p, cr/J) one may predict the qualitative form of its corre- 
sponding phase diagram by using Fig. [5](a). Let us first consider one value of p within 
the range < p < 0.24, where one may have both a critical end point and a tricritical 
point, e.g., the case p = 0.2, discussed previously. Then, for sufficiently low values of cr/J 
one has a phase diagram like the one shown in Fig. [I](b), characterized by a tricritical 
point, a first-order line at low temperatures and a critical end point, associated with a 
phase F2; by increasing a/ J, one reaches the upper limit of region I [represented by the 
continuous line in Fig. [5]^a)] that corresponds to a collapse of the critical end point with 
the zero-temperature axis, like shown in Fig. [Tt^c). Increasing cr/J further, in such a way 
that the probability distribution for the fields still presents a minimum at H = 0, one 
has a tricritical point at finite temperatures [with a typical phase diagram exhibited in 
Fig. [11(d)], and after that, one gets phase diagrams like those shown in Figs. W^e) and (f). 
This sequence of phase diagrams occurs for all values of p within this range. Now, if one 
considers 0.27 < p < p* in Fig. M^a), by increasing a gradually one gets first a typical 
phase diagram characterized by an ordered critical point, like the one shown in Fig. [2](b), 
and then, one may go through the intermediate region between the two dashed lines, with 
a phase diagram as presented in Fig. [5](c), and afterwards, one follows a similar sequence 
of phase diagrams like those shown in Figs. [T](b)-(f). The last two steps of the previous 
sequence of phase diagrams apply if one increases a/ J in the range 0.24 < p < 0.27. By 



19 



varying a j J for p within the range p* < p ^ 0.93 one gets phase diagrams that follow the 
sequence shown in Figs. |3l Finally, for 0.93 ^ p < 1.0 one has a phase diagram typical of 
the Gaussian RFIM, i.e., a continuous phase transition separating the paramagnetic and 
ferromagnetic phases. 

B. Zero- Temperature Critical Frontiers 

Some analytical results may be obtained from the investigation of the free energy and 
magnetization at zero temperature, given by Eqs. ( 12.13^ and ( ]2.14p . respectively. Herein, 
we shall restrict ourselves to a > in Eq. (12. 2p : the case a = 0, i.e., a trimodal PDF 
requires a separate analysis and the results are already known jlJ . fisj ] . 

For that purpose, one applies the same procedure described above for finite tempera- 
tures, starting with the expansion of Eq. (12.141) in powers of m, 

m = aim + asm^ + a^nt' + a^rnJ + O(m^) , (3.7) 



where 
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A continuous critical frontier occurs at zero temperature according to the conditions, 
Oi = 1 and 03 < 0, leading to a relation involving Hq/ J, a/ J and 
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p + (1 — p) exp 



(3.12) 



One notices that, for (Hq/J) = one has (cr/J) = -^Ijix = 0.7979 (Vp), yielding a 
continuous critical frontier for small values of Hq/J. For > 0, one gets a first-order 
critical frontier at zero-temperature, which is usually associated with higher-order critical 
points at finite temperatures. A tricritical point appears at zero-temperature, provided 
that ai = 1, 03 = 0, and < 0, in such a way that. 
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Let us now analyze the coefficient 05 under the conditions ai = 1, 03 = 0. Substituting 
Eqs (Km and (KW in Eq. fIXTU]) . one gets 
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(3.14) 



leading to 05 < for [l — f (^) p] ^ < 3. One should notice that for p = 0, i.e., for a 



double-Gaussian PDF 



the coefficient 05 is always negative; however, in the present 
case one may have 05 = 0, in such a way that one has the possibility of a fourth-order 
critical point. This zero-temperature point is unique for the distribution of Eq. (12.21) and 
it occurs for a set of parameters {p* , Hq/ J,(t* / J), to be determined below. Considering 
as = 0, 



cr 
J 



P 



(3.15) 



and using this result in Eq. (13.13^ . one gets 



H 



° = Vs. 



(3.16) 
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Taking 03 = and using Eqs. fl3.15p and fl3.16l) leads to 

p* = 2 (2 + e^/^^-^ ^ 0.308561 , (3.17) 

which may be substituted in Eqs. fl3.15p . (13.161) to yield, respectively, 

;2 + e3/2)-i ^ 0.369294 ; ^ = 3W- (2 + e^/^)"^ = 0.639637 . (3.18) 

J V TT 



J V27f 



Therefore, for the set of parameters of Eqs. (I3.17P and (I3.18P one has a fourth-order critical 
point at zero temperature, as shown previously in Fig. [2]^d), and also exhibited in the zero- 
temperature phase diagram of Fig.El In the present problem, the zero-temperature fourth- 
order critical point may be interpreted as the threshold for the existence of tricritical 
points in both finite and zero temperatures. For p > p*, there is no tricritical point 
for arbitrary values of a/ J and Hq/J, although an ordered critical point may still occur 
for finite temperatures [cf. Fig. [21(b)], as well as at zero temperature [see Fig. E]. This 
threshold is associated with a very flat PDF, as shown in Fig. JST a), which is in agreement 



with the conditions for the existence of tricritical points UJ, ll2|, Il3 |. 

Zero-temperature phase diagrams of the model are presented in Fig. [6] for several 
values of p. One notices that besides tricritical, and the above-mentioned fourth-order 
point, critical end points and ordered critical points also appear at T = 0. These two later 
points were determined by an analysis of the zero-temperature free energy of Eq. (12.131) . 
similarly to what was done for finite temperatures. Comparing the phase diagrams of 
Fig. [6] with those for finite temperatures, one notices that the width cr produces disorder, 
playing a role at T = that is similar to the temperature; as examples, one sees that 
the phase diagrams for p = 0, p = 0.2, and p = 0.5, in Fig. [6l resemble those shown in 
Figs. W(yd), [D^b), and[3]^b), respectively, if one considers the correspondence a T. 

IV. CONCLUSIONS 



We have investigated a random-field Ising model that consists of an infinite-range- 
interaction Ising ferromagnet in the presence of a random magnetic field following a triple- 
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FIG. 6: Zero-temperature phase diagrams of the infinite-range-interaction ferromagnet in the 
presence of a triple Gaussian random field, for typical values of p. The critical frontiers separate 
the paramagnetic (P) and ferromagnetic (Fi and F2) phases. Critical frontiers and critical points 
are as described in the text. 

Gaussian probability distribution. Such a distribution, which is defined as a superposition 
of three Gaussian distributions with the same width cr, centered at if = and H = 
±Ho, with probabihties p and (1 — p)/2, respectively, is very general and recovers as 
limiting cases, the trimodal, bimodal, and Gaussian probability distributions. We have 
shown that this model is capable of exhibiting a rich variety of multicritical phenomena, 
essentially all types of critical phenomena found in previous RFIM investigations (as far as 
we know). We have obtained several phase diagrams, for finite temperatures, by varying 
the parameters of the corresponding probability distribution, e.g., p and a, with variations 
in p (cr fixed) leading to qualitatively distinct multicritical phenomena, whereas increasing 
cr [p fixed) yields a smearing of an specific type of critical behavior. 

The random-field Ising model defined in terms of a trimodal probability distribution 
[1^ . Q represents, to our knowledge, the previously investigated model that exhibits a 
variety of multicritical phenomena comparable to the one discussed above. Since the 
former represents a particular case (a = 0) of the present model, the parameter a may 
be considered as an additional parameter used herein, when compared with the model 

n n 

of Refs. p^, 115[. It is important to stress the relevance of this additional parameter 
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for the richness of the phase diagrams, as well as for possible physical applications. As 
one illustration of the phase diagrams, one may mention those obtained through the 
zero-temperature analysis of the present model, where one gets a large variety of criti- 
cal frontiers, exhibiting all types of critical points that occur at finite temperatures; this 
should be contrasted with the much simpler zero-temperature phase diagram of the tri- 
modal random-field Ising model (where only an ordered critical point is possible). In 
particular, some of the zero-temperature phase diagrams presented herein resemble those 
for finite temperatures, if one considers the correspondence a ^ T. Therefore, a may 
be identified as a parameter directly related to the disorder in a real system, in such a 
way that an increase in a in the present model should play a similar role to an increase 
in the dilution for a diluted antiferromagnet 26[]. It is important to remind that in the 
identifications of the RFIM with diluted antiferromagnets 0, Q, , the effective random 
field at a given site is expressed always in terms of quantities that vary in both sign and 
magnitude, and so the present distribution for the random fields is more appropriate for 
this purpose. As a direct consequence of this model, the smearing of an specific type of 
critical behavior, due to an increase in the randomness, i.e., an increase in a, or even a 
destruction of a possible multicritical phenomena due to this increase - which represents 
a very common feature in real systems - is essentially reproduced by the present model. 
We have shown that the main effects in the phase diagrams produced by an increase in 
a are: (i) For those phase diagrams presenting two distinct ferromagnetic phases, it re- 
duces the ferromagnetic phase with higher magnetization and destroys the one with lower 
magnetization; (ii) All critical points (tricritical, fourth-order, ordered, and critical end 
points) are pushed towards lower temperatures, leading to a destruction of the first-order 
critical frontiers. 

We have argued that the random-field Ising model defined in terms of a triple-Gaussian 
probability distribution is appropriate for a physical description of some diluted antifer- 
romagnets that in the presence of a uniform external field may exhibit multicritical phe- 
nomena. As an example of a candidate for a physical realization of this model, one has 
the compound FexMgi_xCl2, which has been modeled through an Ising ferromagnet under 
random fields whose distribution appears to be well-represented by a superposition of two 
parts, namely, a trimodal and a continuous contribution [t^. This compound presents a 
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first-order line that disappears due to an increase in the randomness; such an effect has 
been described recently in terms of a simpler model 26|]. However, it is possible that the 



diluted antiferromagnet FexMgi_xCl2, or another similar compound, may present an even 
more complicated critical or multicritical behavior; by adjusting properly the parameters 
of the random-field distribution, the model presented herein should be able to cope with 
such phenomena. 
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